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§1. Introduction 

P ■ 
I— i ■ 

One of the most important characteristic features of a quantum theory, as oppo- 
site to a pure classical theory, is that the set of pure states of a quantum system is not 
a simplex. Consider a quantum system described in a complex Hilbert space TL. Its 
pure states are described rays {e 1 ^ : <j> € [0,2tt]} cH, ||V>|| = 1 or, equivalently, by 
one-dimensional orthogonal projections P = \ip >< Let S be the set of all pure 
states. If [i is a probabilistic measure on S then, according to the standard, linear, 
quantum mechanics, only the density matrix p = J s Pdfj,(P) is observable, not the 
measure \x itself. And yet, as demonstrated in Ref. 1), while searching for the dy- 
namical mechanism of particle track creation and event generation, the measurement 
process introduces a mild non-linearity and leads to a unique random process on S 
- provided the quantum system is coupled, in an appropriate way, to a classical sys- 
tem. **) This mild non-linearity is related to the "time-of-events" observable studied 
fSj ' in Ref. 4) , and calls for a careful re-examination of the fundamental axioms of quan- 

tum mechanism. In particular, if the Event Enhanced Quantum Theory (EEQT) 
developed by Blanchard an the present author will prove to be correct, then precise 
measurements of timing of events may give us more information about the actual 
state of the quantum system than it is usually assumed. The present paper, although 
stemming from the above ideas, deals with a pure formal aspect of EEQT - namely 
with a piecewise deterministic process of jumps (state reductions) on the space of 
pure states in a space with indefinite metric. Spaces with indefinite metric are usu- 
ally being considered as too pathological for applications in quantum theory. 5 ) And 
yet within the formalism of EEQT, where none of the standard probabilistic axioms 
of quantum theory are required, and that because all the interpretation comes from 
a dynamical coupling between a classical and a quantum systems, it is possible to 
have a well defined dynamics with all positive probabilities. In the present paper 
we will follow some of the ideas and notations developed in Ref. 6), though we will 
make an effort at presenting this paper as self-contained. We will start with a linear 
Linblad's type master equation for observables in a Krein space, with a particular 
kind of coupling between a quantum system and a finite-state classical system - cf. 
Eq. (6). We skip the possible unitary Hamiltonian part of the evolution, and we 
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concentrate on the dissipative part alone. By Theorem 2 the coupling generates 
a unique, piecewise deterministic process on the space of pure states of the total 
system. Next, we farther specify the coupling in such a way that we can take the 
trace over the classical system and get the effective piecewise deterministic process 
on the hyperboloid of pure states (of positive norm squared) of the quantum system 
- cf. PDP Algorithm in Section 3. We then examine in details the particular case 
of Krein's space V\ t \ when the space of pure states is isomorphic to the Poincare 
disk. We specify the coupling operators to be fuzzy projections (known also as "spin 
coherent states") as defined in Eq. (25), on a family of symmetrically distributed 
states on the disk (a selection of 33 points from the hyperbolic tiling with Schlafli 
symbol (3,8) - cf. Fig. 1). The plot of the sequence of 100 mln of quantum jumps 
generated by the PDP process suggests that there is a fractal attractor set. Figures 2 
and 3 show the results if the iteration process for two different values of the fuzziness 
parameter e = 0.75, e = 1.75 (for e = 1.0 the projections are sharp). Due to the fact 
that our process is a hybrid one - it has jumps with place-dependent probabilities, 
as in iterated function systems, but it also has dissipative parts of a continues evo- 
lution with a random time length, and also due to the fact that the Poincare disk 
is non-compact, there are no ready to use theorems that we could apply in order to 
ascertain the existence and uniqueness of the invariant measure, as it was done in 
the case of positive definite metric. 7 ) Our numerical simulations suggest a conjecture 
that the probabilities of getting outside a bounded region decrease fast enough with 
the dimension of this region, so that the attractor set and the invariant measure 
supported by this set are well defined. It should be mentioned that, according to the 
linear paradigm of quantum mechanics, mentioned at the beginning of this section, 
the patterns displayed in the figures, although mathematically unique, are physically 
unobservable. According to this paradigm, the only observable is the statistical ma- 
trix resulting from the integration of projections with respect to the measure. In our 
case, due to symmetry of the pattern, it is clear that this statistical state is a mixture 
of the maximal entropy state and the orthogonal projection on the eigenstate of 
(73 corresponding to the eigenvalue 1. It is worthwhile to notice that Lozinski et al. 8 ) 
recently studied iterated function systems on the space of statistical operators of a 
given quantum system. Their approach, even though different in spirit, has several 
areas of overlap with piecewise deterministic processes generated by EEQT. 

§2. Indefinite metric space 

Let V be a complex vector space of finite dimension n with a non degenerate 
sesquilinear form (ip, </>).*^ Let L(V) be the algebra of all linear operators on V 
endowed with the star operation defined by (Aip, <fi) = (ip, A*v). A G L(V) is said to 
be symmetric if A = A*. U € L(V) is said to be unitary if U* = J7 _1 . Let TC(V) 
and U(V) denote the set of all symmetric and all unitary operators resp. Let J{V) 
be the set of all J G H(V) f] U(V) such that the sesquilinear form (ip,4>)j = (tp, J(f>) 



*' It should be noted that the term "Hilbert space" is usually reserved for a space with a positive 
definite scalar product. The term "Krein space" though less known, is more appropriate here. 
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is positive definite. We will assume that J{V) has at least two elements. This is 
equivalent to assuming that the scalar product (ip,<j)) is indeed indefinite, that is 
that there exists a basis {e\, . . . , e^, e^+i, ■ ■ ■ , ek+i} such that {e a , ep) = g a p where 

{g a p} = diag(+l, . . . , +1, -1, . . . , -1) (1) 

with k,l > 1. Whenever we will need to specify the signature of V, we will use the 
notation V^,l) f° r V. There is one-to-one correspondence between elements of J{V) 
and A;-dimensional subspaces of V on which (ip, <p) is positive definite. The group 
Uiy) is isomorphic to U (k, I) and acts transitively on J{V) with the fixpoint group 
isomorphic to U(k) x U(l), thus 

T( V ) = U{kJ) SU{k > l) (2) 

{ ' U(k) x U(l) S{U(k) x 17(0)' 

In particular, for k = I = 2 we get the Cartan symmetric domain A(2) = D4 = 
SO(4, 2)/S (0(4) x 0(2)), with Shilov boundary isomorphic to the compactified 
Minkowski space, while for k = I = 1 we get the Poincare disk A(l) = D\ = 
SO(l,2)/ SO(2) whose Shilov boundary is a circle, with a physical interpretation of 
it being a compactified, circular, time. 

Given J € J{V) we have positive-definite scalar product (ijj, 4>)j = (ijj, J<j>}. Let 
us denote by A J the Hermitian conjugation with respect to this scalar product. We 
call A G C(V) J— Hermitian if A = A J . It is then easy to see that A i-> J A gives a 
one-to-one correspondence between symmetric and J— Hermitian operators. 
Definition 1 We will use the following notation: 

V+ = {^eV: (V,V> >0, (3) 
£ + (V) = {Ae£(V):AV + cV + } (4) 

For V> € V + we will write \\ip\\ for \f (tp, ip). 

Notice that A E C+(V) if and only if PV>|| 2 > for all tp G V+. Each J € J(V) is 
a difference of two orthogonal, complementary projections: 

J = E - F = 2E - I, (5) 

where E = E* = E 2 projects on a /c-dimensional subspace on which {ip, 4>) is positive 
definite, and F = F* = I — E projects on the complementary /-dimensional subspace, 
where {ip, cp) is negative definite. 



§3. Linblad's type semigroups and associated piecewise deterministic 

processes in V(fc,i)- 

It is well known that certain types of quantum dynamical semigroups 9 )' 10 ) are 
in one-to-one correspondence with piecewise deterministic Markov processes on the 
space of pure states of the algebra of observables of the system - cf. 11 )' 12 ) and refer- 
ences therein). By analyzing the proof of the existence and uniqueness theorems it 
is easy to see that these results extend to Vk,i provided the operators that implement 
the jumps are in C + (V). 
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Theorem 2 Let S be a finite set. Let V = Vk,i, and let A be the -k-algebra A = 
® qg5 £(V). A typical element of A is a sequence A = (A a ) ae s, A a G C{V). For 
each pair a, (3 G S, a / (3, let there be given an operator g a p G C + (V). Consider 
semigroup of linear maps 4>t{A) a = exp(Li) : A — ► A, t > 0, with infinitesimal 
generator L defined by 

L{A) a = 9%A P g Pa - ^(A a A a + A a A a ), (6) 

where 

A * = Yl 9p a 9i3a- (7) 

The formula 

tr (PMA) a ) = E / tr(AQ)dp(t; P, a; Q, (3) (8) 

gives then one-to-one correspondence between semigroups <f>t and piecewise deter- 
ministic Markov processes on V + x S with transition probabilities p(t; a, P; (3, dQ) 
described by the following algorithm: 

Suppose that at time to the system is described by a pair (P G P+,ce G S), with P 
bring an orthogonal projection on a vector ip G V+. Choose a uniform random number 
p G [0, 1], and proceed with the continuous time evolution by solving the differential 
equation 

ipt = exp(- -A a ) tp t (9) 
with the initial vector ip until t = t\, where t\ is determined by 

\\M 2 = p (io) 

Then jump. When jumping, change a — > (3 with probability 

Pa~>f3 = ||S73a^iil| 2 /0<i> A a 1p tl ) (H) 

and, if Pa^p ^ 0, change 

V'ti = 9Pa'lpt 1 /\\gpa'lpt 1 \\- 

Repeat the steps substituting t±, ipi,[3 for to, ipo, ct- ft- 

Proof Algebraically, the proof goes, step by step, exactly the same way as in 11 )' 12 ) 
(cf. also 13 ) for the uniqueness), the only thing that we need to check is that positivity 
conditions required in the Hilbert space case also hold when the scalar product is 
indefinite. If f(t) = ||^i|| 2 , then the function f(t) is real analytic, with /(0) = 1 
and f(t) = — (Y>t, A a ^ t ). As long as f(t) > 0, that is as long as ipt £ V + , the deriva- 
tive f(t) is non-positive, because (ipt, A a ip t ) = J2/3^ a \\9a/3^t\\ 2 , and g a/3 G C+(V). 
Therefore / is monotonically decreasing from /(0) = 1 and either it never reaches 
the given value p G (0, 1), so the jump never happens, or, if it reaches the value p 
at a finite time t±, then ip tl G V+ and therefore, because g a p G £+(V), the jump 
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probabilities p a ^p are non-negative. Finally, notice that whenever the denominator 
in Eq. (11) vanishes, this can happen only for discrete values of t\, thus on a set of 
measure zero.D 

Here, as in Ref. 7) we will be interested in a special case of the above process 
that leads to simple iterated function systems with place-dependent probabilities, as 
discussed in. 14 ) Let N he a natural number, and let S = 2 N be the set of all binary 
sequences of length N. Assume C + (V) 3 g a p = gi ^ when a differs from f3 only at 
one, the i-th bit, otherwise g a p = 0. Let 

N 



A = Y.aUi- (12) 



2 = P (14) 



If we take a trace over the index a, we end up with a piecewise deterministic process 
on V+ that can be described as follows: 

PDP Algorithm 1 Start, at time to, with a vector ip € V + . Choose a uniform 
random number p € [0, 1], and proceed with the continuous time evolution by solving 
the differential equation 

4>t = exp(- ±A) Y>t (13) 
with the initial vector ip until t = t\, where t\ is determined by 

Ut 

Then jump. When jumping, change 

lb* -» ipi = gii'tJhii'tA 

with probability 

N 

Repeat the steps replacing to, ipo with t\, ipi. 

§4. Example: - the Poincare disk. 

As a simple example consider the case of k = I = 1. Let us choose Jo € J{V) and 
let ei, e2 be an orthonormal basis in V diagonalizing Jo. Thus (ei,e2) = —^2^2) = 
+1, V can be identified with C 2 - the set of all column vectors 



(16) 



where a\ and 02 are complex numbers, Jq has now the matrix form: 
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while the positive definite scalar product (ip,(j))j is nothing but the standard 

(a,b) Jo = aih + a 2 b 2 (18) 

where the bar in c stands for the complex conjugation of the complex number c. Let 
us denote by the Hermitian conjugate matrix of A. Using our previous notation 
A* = A Ja . The most general Jo — Hermitian matrix, A = A^ is of the form 

A = x°I + <t(x) (19) 

where 

<j(x) = xVi + x 2 a 2 + x 2 a 3 (20) 
and crj are the Pauli matrices: 

i -»—(!; i) < 21 > 

">='•=(?; o) (22) 

"3 = ^ = (n °,1 (23) 



,0, -1, 

and x^ 4 , /i = 0, 1, 2, 3 are real numbers. Note that now Jo is represented by 03. It 
follows that the most general symmetric (with respect to ( , )) operator on V is of 
the form 



cr 3 (x° + <t(x)) = x 3 - ix 2 ai + ix l a 2 + x°cr 3 , (24) 

where x^ are real. In other words moving between symmetric and J— Hermitian 
operators is accomplished by the "Wick rotation" x 1 1— > —ix 2 ,x 2 1— ► ix 1 in the 
parameters x l ,x 2 . 

4.0.1. Fuzzy projections 

In our toy model state of the system is represented by a unit vector^) G V + . 
Proportional vectors describe the same state, therefore it is better to represent states 
by one dimensional orthogonal projections P on positive subspaces. The set of all 
states will be denoted by V+. Quantum jumps will be implemented by operators in 
C + (y). In this paper, as in Ref., 7 ) we will be interested in fuzzy projections of the 
form 

Q(J,e) = \{I + eJ) (25) 

which have the above property. Notice that for e = 1 we have Q(J, 1) = \{I + 
(2Ej — I)) = Ej, and thus Q(J, 1) is a sharp projector. Notice that in V^i we have 
Q(J, 1) € V + because in this case maximal positive subspaces are one dimensional. 
We want to know for which values of e, Q(J,e) is in jC + (V). For this it is enough 
to calculate trace of Q(J, 1)Q(J', e)Q(J, 1), and to find the conditions on e which 
guarantee its positivity for all J € J{V). Let us use the representation using Jo and 
Pauli matrices as above. 
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Definition 3 For all m = (mo, mi, 1712), n = (no,ni,n2) € M 3 set < m • n >= 
mono — mini — m2U2. Let 7i = {m = (mo,mi,m2) :< m • m >= 1, mo > 0}. 
Lemma 4 E'ac/i J G i7(V) is uniquely representable in the form 

J = J(m) = <r(m) = m,o<73 — im20\ + im\U2, (26) 

where m G 7Y. 

Proof J must be symmetric, therefore J = x 3 — ix 2 a\ + ixV 2 + x°<73, with x^ 
real. J 2 = I implies then that x 3 = and (x ) 2 — (x 1 ) 2 — (X2) 2 = 1. Now, 
< ip,Jip >= (ip, JoJip)j , therefore, for < ip,Jip > to be non-negative for all ip, 
Jo J = 03 J = x° + xVi + x 2 o"2 must be positive in the standard sense. That means 
x° > 0. □ 



Definition 5 We will write Q(m,e) for Q(J(m), e), i.e. 

Q(m,e) = i(/ + e<7(m)), (27) 
or, explicitly, in a matrix form: 

0(m, C ) = |f , 1 + l m ° , (28) 
2 \e(— nil + jm2j 1 — emo J 

Note that Tr(Q(m, e)) = 1. 
Theorem 6 (i) For each m, r G H, e G M we have 



where 



and 



Q(m, e)Q(r)Q(m, e) = A(e, m, r)Q(r'), (29) 



w , 1 + e 2 + 2e < m ■ r > 

A(e, m, r) = (30) 



r - 1 , c 2 1 o G H ' ( 31 ) 



/ _ (1 - e 2 )r + 2e(l + e < m • r >)m 
l + e 2 + 2e<m-r> 

(ii) The fuzzy projections Q(m, e) map V + into V + if and only if e > 0. 
Proof (i) Follows by inspection, or by applying the "Wick rotation" to the formulas 
in Ref. 7 ) 

(ii) Let A = A* be a linear operator on V which maps V+ into V+, then we must have 
(Aip, Aip) = (ip,A 2 ip) > for all ip G V with (ip,ip) = 1. Observe that (ip,A 2 ip) = 
Tv(EA 2 E) = Tr(AEA), where E is the orthogonal projection on ip. Since 

A(e, m, r) = Tr(Q(m, e)Q(r)Q(m, e)), (32) 

it is enough to analyze the behavior of the function A. The unitary group of V acts 
on J(V) transitively by its natural action: J ^ UJU*. As the phase of the unitary 
operator cancels out, we can restrict ourselves to the special unitary group, in our 
case SU(1, 1). We have 2 : 1 homomorphism U h-» A(U) of SU(1, 1) onto SO(2, 1) 
given by 

UJ(m)U* = J (Am). (33) 
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The function A in Eq.(30) is SO(2, 1) invariant. Therefore to analyze its behav- 
ior we can always make an 50(2, 1) transformation so that r = (1,0,0), then, as a 
function of m € 7i we get 4A(m, e) = 1 + e 2 + 2emo, which is bounded from below 
only when e > 0. If e > 0, then A(m, e), with m £ H, attains minimum equal to 

at m = (1,0,0). □ 

4.1. Poincare disk model 

The hyperboloid J{V) ~ V+ is isomorphic to the Poincare disk D : 

D = {z e C : \z\ < 1} (34) 

by the isomorphism: 

J(m) = er(m) = mo0"3 + im\a2 — im20\ 44> z (35) 

where m\ + irri2 = 2z/ (1 — \z\ 2 ). 

The action of the unitary group SU (V) on J(V) is isomorphic to the action of 
SU(1, 1) on the Poincare disk D via fractional transformations: A : z i-> (anz + 
012)/ (021^ + 022), the isomorphism being given by ^4 = o^t/o^. 

§5. Piecewise deterministic process on 7-*+. 

As an example consider the piecewise deterministic process described by equa- 
tions (9)-(15), with g, L = <J^Q(m[i], e), i = 1, 2, . . . , N. In, 7 ) where we were studying 

quantum jumps on the complex projective space PC(1) « S 2 , we could neglect com- 
pletely the continuous part of evolution of the state vector, as given by Eq.(9), and 
that because we could choose symmetrically distributed configuration of fuzzy pro- 
jections Q(n[i], e) in such a way that Q(m[i], e) 2 = const • I due to m[i] = 0. 
Such a choice is not possible in our case, with sphere replaced by the positive hy- 
perboloid, because vectors on the positive hyperboloid can not balance to the zero 
vector. Given a sequence m[l], m[2], . . . , m[./V], (m[i], m[i]) = 1, we have 

N 

5^m[i] = Ncm c , (36) 
i=i 

where c is a constant c > and m c , is the centroid vector of unit length square 
(m c ,m c ) = 1. The EEQT algorithm, as described, for instance, in, 11 ) tells us that 
when the quantum jumps are implemented by operators gi then time evolution be- 
tween jumps is implemented by 

^(t) = R(t)4>{0) = e-l^i Q(™M ^(o). (37) 

Due to the fact that <r(m c ) 2 = /, the exponential is easily calculated: 

R(t) = e-K 1+e2 ) [cosh(ect) - sinh(ect)<r(m c )] (38) 
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We will choose our configuration so that m c = (1,0,0), thus <r(m c ) = 03, or, in 
matrix form: 



| (1+e2 ) /e - 



R{t) = e"2 



— ect 

e ect 



(39) 



It is easy to see that when projection Q(r) is the projection on ip(0), Q(r) = 
1-0(0)) (^(0)| then 

\m)\\ 2 = Tr(RWQ(r)), (40) 
which can be easily calculated form the explicit matrix formulas above: 

pit) = \\^{t) || 2 = exp (-t(l + e 2 )) (cosh(2ect) - m sinh(2ect)) . (41) 

For e > and mo > 1 The function 
pit) is monotonically decreasing from 
p(to) = 1 to p(t\) = 0, where to = and 
h = (log ((m + l)/(mo - 1))) /4ec. In 
numerical simulations of the PDP pro- 
cess the time of jump is calculated by se- 
lecting a uniformly distributed random 
number p E (0,1), and then numeri- 
cally solving the equation pit) = p. One 
convenient method of doing it is by in- 
troducing variables y = exp(4ect) and 
s = (1 + e 2 + 2ec)/4ec, and solving nu- 
merically the equation 2py s + y(mo — 
1) — (mo + 1) = for y. uniformly 
distributed random number p G (0,1), 
and then numerically solving the equa- 
tion p(t) = p. One convenient method 
of doing it is by introducing variables 
y = exp(4ect) and s = (l + e 2 + 2ec)/4ec, 
and solving numerically the equation 

2py s + y(mo — 1) — (mo + 1) = for y. To illustrate the resulting iterated function 
system we took N = 33 points m[i] from a hyperbolic tiling with Schlafli symbol 
(8, 3) and plot the results of 100 mln jumps for two e = 0.75 and e = 1.75. 




Fig. 1. 33 points from the hyperbolic tiling 
of the Poincare disk with Schlafli symbol 
(8,3) 
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Fig. 3. Hyperbolic quantum fractal e = 1.75. 
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